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Abstract. Wc address a number of conjectures about the ground state 0(1) loop 
model, computing in particular two infinite series of partial sums of its entries and 
relating them to the enumeration of plane partitions. Our main tool is the use of in- 
tegral formulae for a polynomial solution of the quantum Knizhnik-Zamolodchikov 
equation. 
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1. Introduction 



The present work stems from the investigation of the so-called Razumov-Stroganov 
(RS) conjecture [22] (see also [1, 5]), which is a surprising connection between the 
model of Fully Packed Loops (FPL) and the ground state of the 0(1) loop model. 
Although various corollaries and side results were proved as byproducts of attempts 
to prove the RS conjecture [9, 25], the conjecture itself remains unproven. 

In a series of recent papers [11, 23, 10, 7, 12], it was shown how integral formulae 
for a certain polynomial solution of the quantum Knizhnik-Zamolodchikov equation 
allowed to produce various explicit results, including a connection with the enumer- 
ation of plane partitions (inspired by the conjecture [S]). The same general strategy 
will be used in this article in order to perform some computations on the ground 
state components of 0(1) loop model. We shall address some conjectures of Zuber 
stated in [27]; note that thanks to the Razumov-Stroganov conjecture, they can be 
considered as either conjectures on the FPL model (in which case several of them, 
including one to be discussed below, were proved in [3, 4]), or on the 0(1) loop model, 
the latter point of view being ours. We shall also obtain some new results connecting 
the enumeration of certain classes of Plane Partitions and Non- Intersecting Lattice 
Paths (NILPs) with matchings of the form (7r)p and {pa (see section 2 for an ex- 
planation of the notation), and prove a conjecture of [20]. These can be thought of 
as a small step towards a bijection between Totally Symmetric Self-Complementary 
Plane Partitions (TSSCPPs) and Alternating Sign Matrices (ASMs), the latter being 
in trivial bijection with FPLs, since they provide families of equinumerous classes of 
TSSCPPs and ASMs. 

The paper is organized as follows. In the second section we present the various 
models involved: 0(r) loop model, FPLs, NILPs and TSSCPPs. In the third section 
we describe the quantum Knizhnik-Zamolodchikov equation and the relevant poly- 
nomial solution. In the last two sections we obtain some properties of the entries of 
this polynomial solution indexed by matchings of the form (7r)p and {pa, respectively. 
In particular in each case, we describe the corresponding counting problem for NILPs 



and TSSCPPs. 



2. The models 



In this section we describe the various models that are relevant to this work: 
although the latter is concerned with the 0(1) loop model, we mention here for 
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Figure 1. An example of a state in the 0(r) loop model, here with 
72 = 7. Each row is filled with plaquettes turning left or right. The left 
boundary is identified with the right boundary forming a cylinder. At 
each closed loop is given the weight r. 



motivation the Fully Packed Loop (FPL) model and their connection (the Razumov- 
Stroganov conjecture [22]). We also introduce Non- Intersecting Lattice Paths and 
Plane Partitions. 



2.1. 0[t) loop model. Loop models are an important class of two-dimensional 
statistical lattice models; indeed they present a wide range of critical phenomena 
and many classical models can be mapped into a loop model. Here we consider the 
0(r) loop model. 

Consider a semi-infinite cylinder. Each row is made of 2n plaquettes which can 
contain two possible drawings, as on figure \. We give the weight r at each closed 
loop. 

2.1.1. The space of states. In order to set up a transfer matrix approach, we define 
a state of this model to be the connectivity of the 2n points at the bottom. Figure 
2 represents the connectivity which corresponds to figure 1. The diagrams thus 
obtained are called matchings (or sometimes, link patterns). The number of such 
diagrams of size 2n is the Catalan number c„: 

_ (2n)! 
" n!(n + l)! 

Sometimes it is more convenient to represent states with a series of parentheses. 
The bijection consists in putting a '(' at the start of each arch and a ')' at the end; 
thus, our example becomes: 

()(()()((()))) 
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Figure 2. The same example, but here represented in the form of 
matching (also called link pattern). 




(0,0) (2n,0) 

Figure 3. A state of the 0(r) loop model can be represented by a 
Dyck path. A Dyck path is composed by n NE steps and n SE steps 
making a path above the horizontal line (defined by the first and last 
point). At each '(' corresponds a NE step and at each ')' corresponds 
a SE step. The blue dotted line represent an Young diagram. 

We use the notation (7r)p to represent p parentheses surrounding a matching tt: 

(7r),= (^7r)_^ 
p p 

and ()P for p successive (): 

{Y = {)_^ 

p 

Another way to represent these matchings is using Dyck paths. A Dyck path is a 
path starting at (0, 0) and ending at (2n, 0) composed of n NE steps (or (1, 1)) and n 
SE steps (or (1,-1)), such that the path never goes under the horizontal line defined 
by the extreme points. To construct a Dyck path from a matching we replace each 
opening with a NE step and each closing with a SE step as shown on figure 3. 

Finally, one can also represent a Dyck path as a Young diagram included in the 
(n — 1, n — 2, . . . , 2, 1) Young diagram. The Young diagram is obtained from the Dyck 
path by constructing the complementary of the path. On figure 3 we exemplify how 
to transform a Dyck path into a Young diagram (and vice- versa). 
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qt-q ^Zi 



qt-q ^2 



Figure 4. The two possible plaquettes and their respective probabihties. 




Figure 5. The Yang-Baxter equation. 



In what follows all operators act on the vector space of formal linear combination 
of the matchings |.^) = J^w^-^ 1^)- 

2.1.2. The transfer matrix. The dynamics of the system is encoded in the transfer 
matrix, which represents one row of plaquettes; each plaquette can have a left turn 
or a right turn as explained on figure 4. 

The probability of each plaquette turning left (and in the same way the probability 
of turning right) is defined by a horizontal parameter t and a vertical parameter Zi 
(the i indexes the column). The parameters q and r are related by r = —q — q~^. 
With this parametrization, the weights satisfy the Yang-Baxter equation (figure 5). 
As a consequence, the transfer matrices satisfy the commutation relation 

[T(t\Zi, Z2n),T(t'\zi, Z2n)] = 

SO that, assuming diagonalizability of the T{t\zi, . . . , Z2n), their eigenvectors do not 
depend on the parameter t. We are here especially interested in the ground state 
eigenvector, denoted by 

2.1.3. The Hamiltonian. We are mostly interested in the homogeneous limit where 
all the Zi are equal. In this case one can alternatively define the model using a 
Hamiltonian. 

Definition. The affine Temperley-Lieb algebra. 

Let Ci be the operator that acts as indicated on figure 6, the e2n creates an arch 
between 1 and 2n as our model is defined on a cylinder. This operator obeys the 
affine Temperley-Lieb algebra: 
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Figure 6. The Cj creates a small arch between i and z + 1 as exempli- 
fied in this figure. If the matching tt has already a link between i and 
i + 1 we get the same state but multiplied by a constant ejvr = rvr, as 
the operator Cj leaves vr invariant forming a closed loop. 



CiCj = CjCi if \i — j {mod 2n)\ > 2 



„2 



Consider the Hamiltonian obtained as the logarithmic derivative of the transfer 
matrix T(t|l, . . . , 1), at t = 1, which, up to an additive and a multiplicative constant, 
is given by: 



2n 



(2.1) ^ = -2^e. 

i=l 

2.1.4. The special case r = 1. The case r = 1, or equivalently q = e='=^'^*/^, plays a 
special role. At this value, the Yang-Baxter equation allows us to write [9]: 

(2.2) ( ^5l±l^^i + ^l±l^^e) \^) = \^) 
\qzi-q ^Zi+i qzi-q ^Zi+i J 

where Si permute Zi and Zi+i. Sif{zi,Zi+i) = f{zi+i,Zi). 

Note that in the case Zi = 1 for all i, the ground state of H (and of T) is charac- 
terized by: 

H\tlj) = -2n\^/j) 

Here we use the notation for the specialization at = 1 for all i of |\&). We 
decompose our ground state as = ^7r(^i, • • • , ^2n) k) and similarly for 

2.2. Fully packed loops. A fully packed loop (FPL) consists in a grid n by n, in 
which each point has valence two, i.e. each point is connected to two neighbors. 
We use the boundary conditions defined on figure 7 (equivalent to the domain wall 
boundary conditions of the six- vertex model). 
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Figure 7. A fully packed loop configuration is defined on a grid 
n X n, with paths passing through each vertex. Boundary edges are 
alternatingly occupied or empty; the numbers label the occupied ones. 



If we only consider the pairings between the (labelled) exterior points, we can 
map the configurations of figure 2 to matchings. Note that this is not an injective 
function: generally, there are various FPL states that produce the same matching. 
It is therefore natural to consider the number, denoted A^, of FPL configurations 
with matching vr. The Razumov-Stroganov conjecture (formulated in [22]; see also 
[1, 5]) states the following: 

Conjecture. Let be = '^^'ipn l^r) the ground state of the Hamiltonian 

In 
i=l 

at T = 1, with the normalization condition ipQ^ = 1. Then 

2.3. Non-Intersecting Lattice Paths. These paths are defined on a lattice and 
connect a set of initial points to a set of final points with certain conditions (see Ref. 
[19, 15] for the general framework). The most important feature of Non- Intersecting 
Lattice Paths (NILPs) is that the various paths do not touch one another, this will 
be important on the process of counting them using the Lindstrom-Gessel-Viennot 
(LGV) formula [19, 15]. 

We shall be interested in some classes of NILPs. We pick the area defined by 
X < y < 2x. We define a path as starting at a point [i, 2i) with i a non-negative 
integer less than n. Each path is composed of steps forward and down, ending in 
some point in the line y = x. We are interested in the set of n paths exemplified 
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Figure 8. In the left we see a NILP composed by 5 paths. In the 
middle, we see another NILP, but the p dashed red line are fixed, the 
function that counts these NILPs is called J-'p-i^r+i where n = p + r is 
the number of paths (here 5), p the fixed paths (here 4). In the right 
we find something similar, but the dashed red lines are now vertical 
and we associate them the function Qp^^- 



on figure 8. Each vertical step has a weight r, for example, the first figure in 8 has 
weight r^. 

In the sections 4.5 and 5.2 we will count the number of these NILPs. For this it 
will be important to consider some constraints in the NILPs. We can impose the 
first p paths to be horizontal or to be vertical (we consider the null path at z = 0). 
The associated counting functions will be called Tp-i^r+i and Qp^r, respectively. 



2.4. Totally Symmetric Self— Complementary Plane Partitions. Pictorially, 
a plane partition is a stack of cubes pushed into a corner, gravity pushing them to 
the corner, and drawn in a isometric perspective as exemplified on figure 9. 

Plane partitions were first introduced by MacMahon in 1897. In the pictorial 
representation. Totally Symmetric and Self-Complementary Plane Partitions (TSS- 
CPPs) are Plane Partitions inside a2nx2nx2n cube which are invariant under the 
following symmetries: all permutations of the axes of the cube of size 2n x 2n x 2n; 
and taking the complement, that is putting cubes where they are absent and vice 
versa, and fiipping the resulting set of cubes to form again a Plane Partition. 

As we will see in sections 4.6 and 5.3 there is a bijection between TSSCPPs and 
the class of NILPs described before. More precisely, there is a bijection between 
NILPs with fixed paths and Punctured TSSCPPs. 
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Figure 9. A plane partition. 




Figure 10. The introduced rotation operator p acts in a natural way 
on the space of matchings. 

3. Polynomial solution of the quantum Knizhnik-Zamolodchikov 

equation 

At the value r = 1 of the 0(r) loop model, the ground state eigenvalue takes a 
simple form, and the ground state components become polynomials of the inhomo- 
geneities [9]. We now introduce an equation satisfied by these polynomials. First we 
define the invertible operator p with relations: 



pci = Ci+ip for alH G {1, . . . , 2?T, — 1} 

p2" = l 

In the loop picture this operator corresponds to the rotation operator, see figure 10. 
Consider the following system of equations for homogeneous polynomials \E'^ of 
the variables ^i, . . . , Z2n of degree 5: 
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• The exchange equation (identical to (2.2)): 

(3.1) fl^l±l^^i + -J^±l^e) \^)=sm z = 1, . . . , AT - 1 

• The rotation equation is written in the form: 

p"^^^(Zl, 22, • • • , Z2n) = H'^n{z2, . . . , Z2n, q^ Zi) 

where k is a constant such that = 1, so that K^'^cf^ = 1. 

Although these equations are not what is usually called the quantum Knizhnik- 
Zamolodchikov (introduced in [13] as a g-deformation of the Knizhnik-Zamolodchikov 
equations), it is easily shown that the solution of this system is also solution to the 
(level 1) gKZ equation. 

When q = e^'^'^l'^ (that is r = 1), we obtain the equations that characterize the 
ground state of the 0(1) loop model. In this case the minimal degree for which there 
exists a solution is 5 = nin — 1) (all other polynomial solutions will be multiples of 
this lowest degree solution at g = e^'^*/^). 

In [10, 26], a method is described to construct the solution of this system sys- 
tematically. Schematically, an order is defined on the matchings, and rewrite the 
exchange equation in a triangular form such that it is enough to know for the 
smallest element vr = ()„, in order to be able to compute "^t^^ for all vr. In fact this 
triangular system can be explicitly solved [17, 6]. 

However this method is not particularly convenient for our purposes, and here we 
shall use instead an integral formula for using another basis (presented in section 
3.2). 

3.1. Wheel condition. Using equation (3.1) and knowing that the polynomial has 
degree n{n — 1) we can prove the important relation \I',r(^i, • • • , Z2n)\z-^=q2^.^q^zi ~ 
for all k > j > i, called wheel condition (see [21]). 

The space of homogeneous polynomials in 2n variables of total degree n{n — 1) 
which satisfy the wheel condition has been studied in various articles [16, 21]. The 
interesting fact is that this space has exactly dimension c„ (the Catalan number and 
also the number of different matchings of size 2n), as has been proven in [21] and, in 
a simpler way, in the appendix C of [12]. 

The proof consists in proving that these poljTiomials are completely characterized 
by the values it takes at the c„ points q^ := (g'^\ . . . , q"^"), where = ±1 is another 
way to encode a matching, —1 for the opening '(' and +1 for the closing ')'. In 
what follows, e will represent both the matching and the corresponding series of ±1 
depending on the context. 

In order to proceed we need the following lemma. 
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Lemma 1. In a vector (g*^), we pick one pair = {q~^,q)- If it connects 

+ 1), then we have: 

i~l 2n 

(3.2) vl>^(g^) = - gVO J] " rtO^^l?') 

j=l j=i+2 

where the hat means that we remove the little arch {i,i + 1). Otherwise 

= 

For the proof of this lemma see [26]. This lemma allows us to calculate the value 
of '^■^{q'^) but we shall not need the explicit result. 

3.2. Another basis: the a-basis. Let a = (ai,...,a„) be a strictly increasing 
sequence of integers such that ai > 1 and at < 2i — 1 for all i. We consider the 
following contour integral: 



(3.3) 4.(21,..., 22„) 



l<«<i<2n 



where the integral is performed around the Zi but not around q~'^Zi. 

We claim that these polynomials span our vector space of dimension c„. To prove 
it we need to check several properties: 

• Prove that the $a are indeed, polynomials. It is enough to prove that the 
expression does not have poles, as done in [11]. 

• Check the total degree (5 = n{n — 1)) and the homogeneity. This follows 
directly from the integral formula, knowing that they are polynomials. 

• Check the wheel condition. See [11] for this proof. 

• Prove that the $a are linearly independent. This will be obvious after the 
calculation of the basis transformation between '^t^ and 

3.3. Basis transformation. We assume that the are linearly independent. We 
define the basis transformation as follows: 

(3.4) * (^1, • • • , Z2n) — C*7r,a'^'a(^l, • • • , ^2n) 

a 

And the inverse transformation: 

(3.5) $,(^i,...,Z2n) = Xl^- 

,7r ^ TT 

(^Zl, . . . , Z2n) 
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Note that the polynomials \E'7r can be multiplied by some constant without chang- 
ing their characteristics, so we can define: 

^'0(2:1, 2:2) = $1(21,22) = 1 

3.3.1. Calculating the coefficients Ca,n- We know that these coefficients are unique 
and fully determined by the Cn points {q"). Therefore we only need these Cn points. 
Appendix A of [11] describes how to perform this calculation graphically. We shall 
sketch here the process. For this, we will allow all a = (ai, . . . , an) such that < a^+i 
with 1 < tti < 2i — 1 for all i. 

We pick a little arch in vr, say between i and i + 1. If there are no aj = i the 
is zero. So, a is of the form 

a = (ai, . . • • • ,an) 

k 

After a tedious calculation we obtain: 

i-l 2n 

(3.6) $,(g^) = - q'q^^) H {1 - q-\^^)U,.,^a{q') 

j=l j=i+2 

where a = (ai, . . . , z — 1, . . . , -i — 1, . . . , a„ — 2) and 

k-l 

k —k 

Uk-i = r 

q-q 

is a polynomial of r of degree k. 

Observe that the pre-factor is exactly the same as in (3.2). So, we get: 

The fact that Ci^q = 1 provides an inductive method of calculation. 

There is to our knowledge no direct method to compute C-^^^^a explicitly. We shall 
use an Ansatz later for some entries of C and confirm it by checking it at all values 
of the form (g^). 

3.3.2. Triangularity of the transformation of basis. Now we take some strictly in- 
creasing set a = (ai, . . . , a„) such that < 2i — 1 and the set of openings in the 
matching vr called (tti, . . . , 7r„), in this notation, we see that the two objects are es- 
sentially the same. We define a partial order: vr < a if and only if VTj < for all i. 
We claim that: 

Lemma 2. For Ca^n to be different from zero, ir must be smaller or equal to a. If 
Tii = tti for all i, Ca,n = 1- 
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To prove that we shall use the geometric method presented in section 3.3.1. We 
pick an arch in vr, going between i and j. If we use the geometric method of reduction, 
we see that we need at least (j — i)/2 ak such that i < cik < j , applying this in all 
arches in vr we see that we need n < a. 

If TTj = tti for all i, the calculation of Ca^-K is really simple. 

So, if we set a total order which respects the partial order defined before, Ca^^ is 
a triangular matrix with I's on the diagonal. 

The two facts together prove that the transformation is invertible, therefore the 
$a are linearly independent. 

3.3.3. Coefficients as polynomials of t. The study of these coefficients leads to the 
following property: 

Lemma 3. Let Ya be the Young diagram (Ai, . . . , A,.) corresponding to a ={..., ai = 
Xn-i+i + • • •)? ^'^^ bc thc Young diagram corresponding to the matching tt. 

The coefficients Ca,n are polynomials in t, more precisely: 



where Pa,TT{T) is a polynomial of t with degree baj-K < |^a| — I^ttI — 2. 
We leave the proof of this lemma to appendix C. 

Observe that this lemma remains true for the coefficients C^,a, by using triangu- 
larity of C and the fact that the diagonal elements are one. 



In general, the computation of the polynomials ^.j^ is complicated, and there is no 
general closed formula. But, there are some exceptions. In this section we study the 
polynomials indexed by configurations of the type (vr)p, i.e. a given configuration 
surrounded p parenthesis (see figure 11 for an example). In subsection 4.3 we present 
some properties of the polynomials for high p. 

In all that follows, vr is a link pattern of size 2r, so that (7r)p has size 2n with 
n = r + p. 



(3.7) 




ifY.^^Ya 

lfY. = Ya 

^fY^dY^ 



4. Study of entries of the type (7r)p 



4.1. a-Basis. As done in other articles [11] and explained in 3.3.1 we decompose our 
polynomials in the a-basis. In this case, we find that: 
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Figure 11. A matching (7r)p with p arches surrounding a hnk pattern 

IT. 

Lemma 4. We have the following decomposition 

(4.1) <i/(^^)^{Zi,...,Z2n) = 

,a'^l,...,p,p+ai,...,p+ari.^l: ■ ■ • : ^2n) 

l<ai<...<ar<2r 
ai<2i—l 

where the coefficients are the same that occur in 

(4.2) ^^(zi,...,Z2.)= Yl 

[Zi, . . .,Z2r) 

l<ai<...<ar<2r 
ai <2i—l 

This is, if we find the Cj^^a which satisfy the equation (4.2), these same coefficients 
solve the equation (4.1). 

One checks that the triangularity forces $i,...,p,p+ai,...,p+ar to be written as a linear 
combination of entries of the type "^(-k),,, and these coefficients do not depend on the 
value of p. When one inverts the transformation only these coefficients will matter, 
proving the lemma. 

4.2. Reduction to size r . Now, for a given link pattern tt of size 2r , we can calculate 
polynomials ^/'(7r)p for all p. For example: 



T 



(4.3) ^(000)^ = ^r(2r2n2 - bnr^ + Sr^ + 6) 

(n-2)(n-l) 2 
^((())())p = 2 ^ 

where r = —(q + q~^)- These are the steps of the calculation: 

• For given p and r we compute the Ca,n for all vr and a of size r. 

• We invert this relation. 
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• We calculate 



^l,...,p,p+ai 



,...,p+aAzi, . . . , Z2n) = (-l)^^) JJ {qZi - q ^Zj) Yl ((l^i ^ ^ ^^j) 



l<«<j<p p<i<j<2n 




rT dwiUi<i<j<ri^j-Wi)iqWi-q ^Wj)Yli<j<p{(lZj - q V) 



(4.4) 




'l,...,p,p+ai,...,p+ar 




1 + TUm + M«Mm)(l + 



4.3. Expansion for high p. Zuber conjectured the polynomial dependence in p and 
large p behavior of the number of Fully Packed Loop configurations with connectivity 
(7r)p (conjecture 6 in [27]). This was subsequently proved in [4]. Alternatively, due to 
the Razumov-Stroganov conjecture, one can expect the same behavior for the ground 
state entries of the 0(1) loop model. Here we generalize it to the gKZ solution for 
any r: 

Theorem. For matchings of the type {tt)p, the polynomials ipT^ can be written in the 
following form: 



where Y is the Young diagram defined by ir, \Y\ its number of boxes, and PyiT^n) is 
a polynomial in n and r of degree \Y\ in each variable with integer coefficients. 
In the limit of high n, the polynomials behave like 



where dimy is the dimension of the irreducible representation of the symmetric group 
associated to Y . 

As the basis transformation is triangular, we can write: 



where a is equivalent to vr and h < a means that the Young diagram of h is inside of 
the Young diagram of a. We denote the Young diagram corresponding to a by Ya- 






b<a 
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Note that, by lemma 3, C,r,6 are polynomials of r with integer coefficients and 
exponent no more than \Ya\ — \ Yb\ — 2. 

We now prove that the integral of the ffist term, corresponding to the largest Young 
diagram in the decomposition, is a polynomial of n with the asymptotic behavior of 
the theorem. The other terms will possess the same polynomiality property, and they 
will be of lower degree in n and r (noting that the power of r in the non-diagonal 
elements is less than ly^l — \Yb\ and does not affect our conclusion). 

We want to calculate the integral 



/f dill 

i=l * i<j<r 



We replace the term ni<j(l + ™i + ^i^j) with ni<i(-'- + ^^i) = 11^(1 + 



TU, 



because any term with UiUj in the product is formally identical to the contribution 
of a smaller diagram. We then compute 

n2l!i^<^+™-'"""n(%-".) 

i=l * j>i 

-i:(-ir/-/n^(i + ™.r-'<-' 

"■iE(-i)°ri^""''^'"' 



i=l 



where Xn-i+i = Oj — i is the size of each row in the Young diagram and (—1)'^ is 
the sign of the permutation a. We obtain that the coefficients can be written as a 
sum of integers divided by Ylii'^ ~^ ■^i ~ '^iV- which divide as a consequence of 

E.(^ + A,-a,) = |y|. 

The dominant contribution ClS db function of n is 



r 



^1 



i+A„-i+i— (Ti J" 1 



A (z + A„_i+i - cTj)! ^ fJi + Ki-i+i - cri)\ 

dimF 



\Y\ 



where the last equality can be found, for example, in the fourth chapter of [14]. 
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4.4. Sum rule. Pick an a of the form (ai, . . . , a,.), with = 2z — 1 or = 2'i — 2.^ 
Let £(a) be the set of matchings whose openings on odd sites are exactly the odd 
elements in a. 

From [11], section 3.3, we know: 



7rG£(a) 



It follows that, using lemma 4: 



(4.6) r^'-^»'^'$a = $^r°-vl/(^)^ 

a=(l,...,p,p+ai,...,p+ar) ^ 
ai=2i—l or 2i— 2 

where o^r counts the openings in even sites, of the matching vr, and (r^ — ^ • a^) 
counts the number of even in a, as o^r = — if vr G £(a). 

In (4.6), at = 1 for all z, the l.h.s. is equal to: 

• • • r n 2^^^2f-i n ~ "'^(■^ + + M/tim,)(i + rM„)p+^ 

1=1 I l<rn<r 

but this is exactly 

p+i 

We can now state the main result of this section: 

Theorem. Let count the number of arches of vr opening at an even site, we have 
the following result: 

(4-7) Yl = 

IT of size 2r 

At r = 1, we get the proof of conjecture 8.i of [27], but for the 0(1) loop model: 
(4-8) = AiYUi 

n of size 2r 

At p = 0, we can use the rotation symmetry and obtain the already known formula: 

TT of size 2r 



Observe that if oi = 0, = 0. 
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r 




Figure 12. In order to apply the LGV formula we label the starting 
points of the paths by 1, . . . , r and the ending points by Li = 1, . . . , L^. 
As the paths do not intersect each others Lj+i > Lj. 



4.5. A NILP formula. We can interpret the result of the previous subsection in 
terms of the NILPs previously defined, cf figure 8. We fix p paths as shown in 12, 
and using the LGV formula, we are able to calculate the number of different NILPs. 
We label the paths with i = 1, . . . , r and call the final locations Li. 
If we only consider one path going from i to Li. 




we require that Li = 1. We give a weight r to each vertical step. 
The LGV formula tell us that the number of paths is equal to 




We can use a contour integral form: 
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SO, the number of NILPs can expressed by: 
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p,r 



TT^ y det 

i=\ \=L\<...<Lr 



|P+fe— 1 



2/c — Ij^ 



1 = 1 * l=Ll<...<Lr 

i=l * i<j 



Uj — Ui 



l<fc,j<r 



l<fc,jr<r 



UiUj 



n 



where the equahty between the second and the third hne can be found, for example, 
in the fourth chapter of [2]. 

Applying the formula obtained in appendix D of [12], we transform the formula 
into: 

(4.9) J^p,, = / • • • / n 2^ki^'-^ ^ 'Tu.Yi^ + Ui) - + ™i + UiUj) 

i=l * i<j 

We recognize this integral formula: 

TT of size 2r 

There does not seem to be a simple closed formula for jFp even at r = 1, as 
suggested by the expressions of [27] for small values of p. 



4.6. Punctured— TSSCPPs. Recall that there is a bijection between NILPs and 
TSSCPPs. As is seen on figure 13, fixing the first p+1 paths to be horizontal 
amounts to fixing the central hexagon of size 2p. Observe in addition that the 
triangle which on the figure contains the paths is a fundamental domain i.e. defines 
the whole TSSCPP. Consequently, Tp^r also counts the number of TSSCPPs with 
weight r for each blue face (the faces containing a vertical step of the paths) in the 
fundamental domain. See [7] for a similar interpretation of partial sums in a related 
model in terms of punctured plane partitions. 
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Figure 13. This TSSCPP with fixed central hexagon, of size 2p, cor- 
responds to the NILPs of figure 12 with p + 1 fixed horizontal paths. 
The corresponding NILPs are drawn in the fundamental domain (the 
triangle between the two blue dotted lines). 

5. Study of entries of the type (^a 

In this section we consider entries ■j/'tt of the polynomial solution of gKZ (in the 
homogeneous limit) which correspond to the matchings with p openings at the be- 
ginning. Explicitly, such a matching is of the form (^a, where a is not a state but 
a sequence which contains n closings and r = n — p openings, see figure 14 for an 
example. Equivalently these are the matchings for which the p first points are not 
connected to each other. 

5.1. Basis transformation. We are here interested in the sum of such polynomials 
a- Recall that there is no systematic way to obtain the corresponding expression 
in the a-basis. However, based on some numerical experimentation, the following 
formula can be guessed: 
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p a 

Figure 14. A state {pa is a state that has at least p openings at the 
left. 

Lemma 5. The sum of the polynomial of the type "^{pa is: 

(5-1) 5^*(p-= $1.2,...,p,a„...,a. 

p<ai<2p+2i~l 

for all ai even 

The proof consists in evaluating this equality at all c„ points of the type (g^). 
It is obvious that the l.h.s. evaluated at (g^) such that e can not be written as {pa 
is zero. 

As to the r.h.s., we integrate the first p variables: 
^ ( 1^m TT ^ -1 ^ / / TT ]\J>^{wj - w,){qwi - q-lw^) 

=(-!)<'' n - ''"''^) n fe^' - ''"''^) 

i<j'^^P p<i<j<2n 

dwi - - q~^Wj) ]\^^i{qzi - q'^Wj) 



X 



n 



where . ■ . ,bp, bp+i, . . . , 6„} = {1, . . . ,p, Oi, . . . , a„}. This sum is zero for all {q") 
such that e do not have p openings at the left. 

We now proceed by induction on r. If r = n — j» = we find that both sides are: 

(5.2) 

We now want to show that for all e of the type {pa the r.h.s. satisfies the same 
recurrence as the l.h.s. Take an e which has a pairing ("little arch") {i,i + 1). Using 
lemma 1, we rewrite the l.h.s. as: 

i-l 2n 

(5.3) Yl ^u^') = ^"^""'^ U(^ ' - ^'^^^) n - ^"'^^^) E "^Uq") 
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where the hat means that we remove the httle arch from (^a^ and e. It is obvious 
that (pd passes by all matchings with n — 1 arches and p openings at the beginning. 

We proceed similarly with the r.h.s. We pick the same vector [q'^). We suppose 
that i is even (for i odd the reasoning is the same). 

If z ^ {dj}, the expression vanishes. So we pick aj = i, and as i is even we have 
ttj+i > i + 1. When we integrate on Wj in expression (3.3), by the rules defined 
in 3.3.1, we obtain a similar formula with a reduced vector e of size n — 1 that is 
obtained by removal of the little arch. 

The integral formula is modified in the following manner: 



i-l 2n 

ai+i>ai j=l j=i+2 
p<ai<2p+2i-l 

Oi+i^ai+l for allai even \ ^ ^ / e\ 

X 2^ ^l,2,...,pM,-,&r-Aq ) 

a.i+l>a,i 
p<di<2p+2i-l 
Si+l^ai+lfor all di even 



where 



(5.4) di 



di i < j 

flj+i — 2 otherwise 



This proves the lemma. □ 
Now we can calculate the limit Zi = 1 for all i. Using the change of variables 



Wi-1 
Ui = 



qwi-q ^ 



and integrating the first p variables, we obtain: 



a ai+i>ai m=l "™ 

p<ai<2i— 1 

a^+i^ai+lfor alia, even JJ _ + + 

1<'<'TI<'' 



if a does not have a little arch (i, i + 1), the term is zero. 
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To sum over all possible a, we can consider only the odd Oj and multiply by (l + n^), 
simplifying in this way the conditions: 



^ . , , . ^ J ^ — 1 



dUm {I + TUmYil + Um) 



p<ai<2i-l 

^iodd II (^Um-Ui){l + TUm + UmUl) 

l<l<m<r 

We write bi = p — (ai + l)/2 + i: 



a 0<b,+i<fei m=l 



(iw^n (1 + TUmY{l + Um) 2b„ 

2m u^rt'""-' """^ 



Y\ {Um -Ul){l + rUm + UmUi) 
l<l<m<r 

The upper bound on bi was relaxed, since it only excludes zero terms. 
A standard calculation gives us the formula: 

1 



2bi 



Replacing, we get: 



,r- r-^ Sr^ , / 1 Tl ^^rn (1 + + Mm) 

J]^ {Um -Ui){l + TUm + Mm^/) 

l<.l<m<.r 

5.2. NILPs. Instead of fixing the first p paths to be horizontal, as in section 4.5, 
we fix them to be vertical. We count the number of paths with a weight r for each 
vertical step (the first p fixed paths excluded). 

We use the same method that was used for the '^{■k)p polynomials. On figure 15 
we show how to label the entries. 

Consider all the paths going from i to Lj, we get: 



' I, Li — ' 

Or, in a contour integral form: 



p -\- i — 1 
p + 2i- Li-1 



_ f du {1 + Tu)P+'-^ 
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Figure 15. To apply the LGV formula we need to label the starting 
points of the paths by 1, . . . , r and the ending points by 1 < Li, . . . , L^. 
As they do not intersect each others Lj+i > Lj. 



We apply the LGV formula: 



27ii ^ ^ 

i=l l=Li<...<Lr 



(1 + rMfc)f+'^-^ 



p+2k~Li 



U, 



n 



dUi (1 + TUi) 

2m 



E det 



l<fc j<r 



-,-r du, (1 + TU^)P+'-' 1 -,-r , 

-l--l-27rz 111 



— Ui 



1=1 



n 

i=l 



2m 



u 



p+2i-l 



We now use the following identity (similar to the one formulated in [11] and proved 
in [24]), proved in appendix B: 

fr du, Ui<jiuj -Ui){l + TUj + UiUj) 

n ^ .r-(i-nw^i) ■ ■ ■ ■ ^ 

where . . . , Ur) is some symmetric function, here fl{ui, . . . , Ur) = 111=1(1 + 

TUiY (l+Ui) , without poles in the integration region. We thus obtain exactly equation 
(5.5). 
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In [18], Krattenthaler gave an explicit formula for Qp^r at r = 1: 

n::; (3.a::SL.. Uto'^^^^p + 2. + m2^)l a r is even 



(5.6) Qp^r 



(3p+2j+l)!{p+2j)! lli=0 

2^ Kzi (3,5:;;s^2.)i ut\''^\2p + 2^)1(2. - 1)! if r is odd 



Remarkably, these formulae coincide with those conjectured in [20] for J2a "^(p" at r = 
1 (more precisely, what was conjectured, cf their Eqs. (40-42), was the probability 
that p consecutive points are disconnected from each other in the 0(1) loop model, 
that is the ratio of V'(pa) by the full sum). Correcting a misprint in their Eq. (42) 
we have 

(5.7) g,.^^(2(, + .),,) 



(5.8) S{L,p) 



S{2p,p) 



p even 
p odd 



The equality of (5.6) and (5.7) can be obtained by direct computation, treating 
separately the parities of p and of r. 



5.3. Punctured— TSSCPPs. As is seen on figure 16, fixing the first p paths amounts 
to fixing a central hexagonal star of size p. Consequently, Qp^r also counts the num- 
ber of TSSCPPs with weight r for each blue face in the fundamental domain (the 
triangle outside the hexagonal star, which on the figure contains the paths). 

Appendix A. Numerical data 

In this section we give some results, that confirm the theorem in 4.3, showing 
explicitly the form of the polynomials. We use r = 1 for simplicity.^ 
For r = 2: 

^(0{))p =in - 1) 

For r = 3: 

(n - 2) , , 

{n-2)in-l) 
^((0)0)p = 2 



'see (4.3) for some examples with general r. 
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Figure 16. This TSSCPP with fixed central hexagonal star, of size 
p, corresponds to the NILP in figure 15 with p fixed vertical paths. 



r = 4: 



(OOOOp 

(((()))())p 
({()()){))p 
({()){()))p 
((())()())p 
({){()){))p 



' ^ {4n^ - 32n^ + ISSri^ - 394n + 540) 



180 

n - 3)(ra - 2){n - 1) 



6 

n — 3)(n — 2) 2 



24 

n-3)(n-2)2(n- 1) 



(3n^ - 7n + 16) 



12 

n - 3){n - 2){n - 1), 2 



24 



n - 4n + 12) 



n — 3) 

^ ^n^ - Tn^ + 24n2 - 48n + 60) 
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Appendix B. Proof of an antisymmetrization identity 

In this section we will follow the procedure of Zeilberger [2 1] . We want to prove 
the following equality: 



UjUi) 



where A is the antisymmetrization operation on the variables ui, . . . ,Ur, and the 
subscript < means that we only are considering the monomials of the kind oc Yl '^T 
with flj < for all i. 

We can directly antisymmetrize the right term: 

where A(Mj~^) = nj->i(^7^ ~ ^j"^) ^^e Vandermonde determinant. 

In what follows we shall call the two sides of equality (B.l) before the truncation 
with <, Ap r{ui, . . . , u,.) and Bp r{ui, . . . , u,.) respectively. 

The proof will be done by induction. The first step is to calculate the case r = 1: 



A 



< 



Next, we suppose that Ap^r~i{ui, . . . , Ur- 1)< = -Bp,r-i(^i, • • • , Ur^i)<. We have 



Using the hypothesis and the fact that < is a linear operator, we rewrite the 
conjecture as: 

J \ Bp^r-liUi, ...,Uj,...,Ur 

Working a little bit the expression, we obtain: 

This equality is a consequence of the following identity, which was pointed out in 
[24]: 

11^(1 ~ UjUj) TT 1 + TUj + UjUj _ 
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In order to prove this identity, we replace uj — > uj^ for all j: 
Or, written in another form: 



I L-^,.-(u,- - Ui) I |,(r + n,; + Ui) v . v 

We recall how to prove this identity using the Lagrange interpolation formula: 



Theorem. Let he P{z) a N—1 (or less) degree polynomial in z and let he (wi, . . . , wjv^ 
different points. So these points define the polynomial that can he written hy: 

Z - Wi 



p{z) = J2p{w,)11 



. Wj - Wi 



Corollary. The maximal coefficient of P{z) is: 

j2 ^(^^•) 



Let be a and (3 the two roots of (1 + tuj + u'^)- Let P be 

P{z) = (r + 2z) + Ui{T + ^))(1 - u,z) 

i 

It is a polynomial of degree (2r + 1), so by the Lagrange interpolation formula, and 
using the points . . . ,Ur, —ui — r, . . . , —Ur — r, a, [3) to describe P{z), we obtain 
the formula: 

2(-iY- rr = V ^^"^^'^ 



P(-r - 



r — M,- — a 

3 ' 



+ 
+ 



)(-r - - /?) ]\i{-T -Ui- Uj) ]\i^^{-Uj -T + Ui + r) 
P(a) 



[a - (3) Yiiia - Ui){a + + r) 



{^3-a)U^{P-u.){(3 + u, + T) 

Using P{—T — Uj) = —P{uj), a + P = —r and aP = 1, we observe that the first 
two terms are identical (and identical to the l.h.s. of (B.3)), while the sum of the 
last two terms simplifies to 2(— l)*". Thus, we get (B.3). Proving this equality. 
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Figure 17. We represent the Dyck paths corresponding to a = 
(1,3,5,6,7) (the black hne, under) and tt = ((()(()))) (the red hne, 
above). The openings of tt are marked with blue little circles and the 
closings with green little squares. 

We can rewrite the main equality as a contour integral formula: 

ndui (1 + rui)'-^ ]\j>^{uj - Ui) 

where f2(Mi, . . . , m^) is a symmetric function in all Ui without poles in the integration 
region around Mj = 0. 

Appendix C. Proof of lemma 3 

In this section we intend to sketch the proof of lemma 3 about the value of C^^^. 

To each a and tt we associate a Young diagram Ya and Ftt, if Ya does not contain 
^TT, Ca,-K = as explained in section 3.3.2. If Ya = Y^, we obtain, trivially, Ca^^ = 1. 
The interesting case is when C Ya. 

For illustration purposes we shall use an example: let a = (1,3,5,6,7) and 
TT = ((()(()))), with the associated Young diagrams Ya = (2,2,2,1,0) and Y^, = 
(1, 1, 0, 0, 0). We can also represent them in the form of Dyck paths as on figure 17. 

As shown in section 3.3.1, each link contributes with a factor Uk- If a link starts 
at r and finishes at s, k is given by: 

S — T" -j- 1 

k = tl{aj such that r < < s} 

Or, in a Dyck path representation, the k is given by a counting in the NE and SE 
steps in the path corresponding to a, more precisely: 

jjlNE steps between r and s} — MSE steps between r and s} — 1 
k = 
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which is the same as counting squares of the skew Young Diagram Yaj-n'- those under 
the opening with sign plus and those under the closing with sign minus (as on figure 
17). In the example we get: f/if/g ~ where the \J\ corresponds to the arch 
between positions 5 and 8. 

If we ignore the individual arches, we note that the maximal exponent is precisely 
the sum all the squares (with the sign ±), in the example this number is 3 — 2 = 1. 

Knowing that there is always at least one square with a minus sign, we verify that 
the maximal exponent of r is lYal — I^ttI — 2. 
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